home *** CD-ROM | disk | FTP | other *** search
/ Sprite 1984 - 1993 / Sprite 1984 - 1993.iso / src / lib / c / etc / RCS / isinf.c,v < prev    next >
Text File  |  1992-03-27  |  3KB  |  177 lines
  1. head     1.2;
  2. branch   ;
  3. access   ;
  4. symbols  sprited:1.1.1;
  5. locks    ; strict;
  6. comment  @ * @;
  7.  
  8.  
  9. 1.2
  10. date     92.03.27.13.37.14;  author rab;  state Exp;
  11. branches ;
  12. next     1.1;
  13.  
  14. 1.1
  15. date     90.11.02.07.44.01;  author rab;  state Exp;
  16. branches 1.1.1.1;
  17. next     ;
  18.  
  19. 1.1.1.1
  20. date     91.12.02.21.43.17;  author kupfer;  state Exp;
  21. branches ;
  22. next     ;
  23.  
  24.  
  25. desc
  26. @@
  27.  
  28.  
  29. 1.2
  30. log
  31. @Fixed some lint errors.
  32. @
  33. text
  34. @/* 
  35.  * isinf.c --
  36.  *
  37.  *    Machine-dependent procedure to determine whether a double is 
  38.  *    infinity.
  39.  *
  40.  * Copyright 1989 Regents of the University of California
  41.  * Permission to use, copy, modify, and distribute this
  42.  * software and its documentation for any purpose and without
  43.  * fee is hereby granted, provided that the above copyright
  44.  * notice appear in all copies.  The University of California
  45.  * makes no representations about the suitability of this
  46.  * software for any purpose.  It is provided "as is" without
  47.  * express or implied warranty.
  48.  */
  49.  
  50. #ifndef lint
  51. static char rcsid[] = "$Header: /sprite/src/lib/c/etc/RCS/isinf.c,v 1.1 90/11/02 07:44:01 rab Exp Locker: rab $ SPRITE (Berkeley)";
  52. #endif /* not lint */
  53.  
  54. #include <math.h>
  55. #include <machparam.h>
  56.  
  57. #if BYTE_ORDER==BIG_ENDIAN     
  58. #define MSW 0
  59. #define LSW 1
  60. #endif
  61. #if BYTE_ORDER==LITTLE_ENDIAN
  62. #define MSW 1
  63. #define LSW 0
  64. #endif
  65.  
  66.  
  67. /*
  68.  *----------------------------------------------------------------------
  69.  *
  70.  * isnan --
  71.  *
  72.  *    Return whether a double is equivalent to infinity.
  73.  *
  74.  * Results:
  75.  *    1 if the number is infinity, 0 otherwise.
  76.  *
  77.  * Side effects:
  78.  *    None.
  79.  *
  80.  *----------------------------------------------------------------------
  81.  */
  82.  
  83. int
  84. isinf(value)
  85.     double value;
  86. {
  87.     union {
  88.     double d;
  89.     long l[2];
  90.     } u;
  91.  
  92.     /*
  93.      * Put the value into a union so we can check out the bits.
  94.      */
  95.     u.d = value;
  96.  
  97.  
  98.     /*
  99.      * An IEEE Std 754 double precision floating point number
  100.      * has the following format:
  101.      *
  102.      *      1  bit       -- sign of Mantissa
  103.      *      11 bits      -- exponent
  104.      *      52 bits      -- Mantissa
  105.      *
  106.      * If the exponent has all bits set, the value is not a 
  107.      * real number.
  108.      *
  109.      * If the Mantissa is zero then the value is infinity, which
  110.      * is the result of division by zero, or overflow.
  111.      *
  112.      * If the Mantissa is non-zero the value is not a number (NaN).
  113.      * NaN can be generated by dividing zero by itself, taking the
  114.      * logarithm of a negative number, etc.
  115.      */
  116.  
  117.     /*
  118.      * check the exponent
  119.      */
  120.     if ((u.l[MSW] & 0x7ff00000) == 0x7ff00000) {
  121.     /*
  122.      * See if the Mantissa is zero.
  123.      */
  124.     if ((u.l[MSW] & ~0xfff00000) == 0 && u.l[LSW] == 0) {
  125.         /*
  126.          * Infinity.
  127.          */
  128.         return (1);
  129.     } else {
  130.         /*
  131.          * NaN.
  132.          */
  133.         return(0);
  134.     }
  135.     } else {
  136.     /*
  137.      * Normal.
  138.      */
  139.     return (0);
  140.     }
  141. }
  142.  
  143. @
  144.  
  145.  
  146. 1.1
  147. log
  148. @Initial revision
  149. @
  150. text
  151. @d18 1
  152. a18 1
  153. static char rcsid[] = "$Header: /sprite/src/lib/c/etc/sun3.md/RCS/isinf.c,v 1.2 90/09/11 14:47:35 kupfer Exp $ SPRITE (Berkeley)";
  154. d22 1
  155. d24 9
  156. d87 1
  157. a87 2
  158.     if ((u.l[0] & 0x7ff00000) == 0x7ff00000) {
  159.  
  160. d91 1
  161. a91 1
  162.     if ((u.l[0] & ~0xfff00000) == 0 && u.l[1] == 0) {
  163. a108 1
  164.  
  165. @
  166.  
  167.  
  168. 1.1.1.1
  169. log
  170. @Initial branch for Sprite server.
  171. @
  172. text
  173. @d18 1
  174. a18 1
  175. static char rcsid[] = "$Header: /sprite/src/lib/c/etc/RCS/isinf.c,v 1.1 90/11/02 07:44:01 rab Exp Locker: rab $ SPRITE (Berkeley)";
  176. @
  177.